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ABSTRACT 

The most X-ray luminous cluster known, RXJ1347-1145 (z = 0.45), has been the object of extensive 
study across the electromagnetic spectrum. We have imaged the Sunyaev-Zel'dovich Effect (SZE) at 90 
GHz (A = 3.3 mm) in RXJ1347-1145 at 10" resolution with the 64-pixel MUSTANG bolometer array 
on the Green Bank Telescope (GBT), confirming a previously reported strong, localized enhancement 
of the SZE 20" to the South-East of the center of X-ray emission. This enhancement of the SZE has 
been interpreted as shock- heated (> 20keV) gas caused by an ongoing major (low mass-ratio) merger 
event. Our data support this interpretation. We also detect a pronounced asymmetry in the projected 
cluster pressure profile, with the pressure just east of the cluster core ~ 1.6 x higher than just to the 
west. This is the highest resolution image of the SZE made to date. 

Subject headings: galaxies: clusters: individual (RXJ1347-1145); cosmology: observations; cosmic 
microwave background; GBT 



1. INTRODUCTION 



The rich cluster RXJ1347-1145 (z = 0.45) is the most 
X-ray luminous galaxy clust er known (Schindler et al. 
1995] |1997| |Allen et al J2002J ) and has been the object of 



extensive study at radio, mil limeter, submillimeter , op- 
tical and X-ray wavelengths ([K itayam a et al. 2004[ Ko- 
' matsu et al.|[200l| [Gitti et al.||^0(j71 |Allen et al.|[ 2002 ; 



Schindler""etaTni997; 



et al.||2008[ ICohen 



Pointecoutea u et al.||1999 |Ota| 
Kneib||2002| |Bradac et al.| |2008 ; 
Miranda et al |I2QQ8| ). discovered in the ROSAT All- 
Sky Survey, RXJ 1347-1145 was origin ally thought to 
be a dyna mically old, relaxed system ( |Schindler et al.| 
1995 1997[ ) based on its smooth, strongly- peaked X-ray 
morphology — a prototypical relaxed "cooling- flow" clus- 
ter. The NOBA 7 b olometer system on the 45-meter 



Nobeyama telescope (Kitayama et al. 2004 Komatsu 



et al. 2001 ) has made high-resolution observations (13" 
FWHM, smoothed to ~ 19" in the presented map) of the 
Sunyaev-Zel'dovich effect (SZE) at 150 GHz which indi- 
cate a strong enhancement of the SZ effect 20" (170 kpc) 
to the south-east of the peak of the X-ray emission, how- 
ever. Hints of this asymmetry had been seen in earlier, 
lower resolution measurements with the Diabolo 2.1 mm 



photo meter on the IRAM 30-m (Pointecouteau et al. 
1999). The enhancement has been interpreted as be- 
mg due to hot (T e > 20keV) gas which is more dif- 
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ficult to detect using X-rays than cooler gas is, owing 
to the lower responsivities of imaging X-ray telescopes 
such as Chandra and XMM at energies above ~ 10 keV. 
In contrast, the SZE intensity is proportional to T e up 
to arbitrarily high temperatures, aside from relativistic 
corrections which are weak at 90 GHz, so such hot gas 
stands out. The feature is consistent with the presence 
of a large substructure of gas in the intra-cluster medium 
(ICM) shock-heated by a merger, as is seen in the "Bullet 
Cluster" 1E0657-56 flMarkevitch et al.||2002) ; this inter- 



pret aion has been sup ported b y more recent observations 



(e.g. |Allen et al. 2002||Ota et al.|2008[ ). Thus, rather than 
being an example of a hydrostatic, relaxed system, high- 
resolution SZE observations suggest that the observed 
properties of the ICM in RXJ 1347- 1145 are strongly af- 
fected by an ongoing merger. This is a striking caution- 
ary t ale for ongoing blind SZE surveys (|Carlstrom et al.| 



2002), for which useful X-ray data will be ditticult or im- 



possible to obtain for many high-z systems, as well as 
a sign that our current understanding of nearby, well- 
studied X-r ay clusters may be dramatically incomplete. 

Reports (|Komatsu et aLl|2001| |Pointecouteau et al. 
200lf |Kitay ama et al.l 120041) of a strong enhancement 



of the SZE away from the cluster center are based on 
relatively low-resolution images compared to the size of 
the offsets and features involved. SZE images at lower 
frequencies also show show substantial offsets between 
the p eak of X-ray and SZE e mission; for instance, the 21 
GHz flKomatsu et al.||2001[ ) and SZE peak is - 20" to 



the S E ol the X-ray peak, and the 30 GHz ( |Reese et al. 
2002[ ) SZE peak is - 13" to the SE of the X-ray peak. 
The situation is further complicated by the presence of a 
radio source in the center of the cluster. We have sought 
to test these claims, and to begin to untangle the astro- 
physics of this interesting system, with higher resolution 
imaging at a complementary frequency. In this paper 
we present the highest angular resolution image of the 
SZE yet made. We observed RXJ1347-1145 with the 
MUSTANG 90 GHz bolometer array on the Robert C. 
Byrd Green Bank Telescope (GBT). At the redshift of 



2 



the cluster (and assuming Q\ = 0.3, Qtot = l,h = 0.73) 
the GBT+MUSTANG 9" beam corresponds to a pro- 
jected length of 54 kpc. The observations are described 
in § [2] and the data reduction in § [31 Our interpretation 
and conclusions are presented in §^4] 

2. INSTRUMENT & OBSERVATIONS 

MUSTANG is a 64 pixel T ES bolometer array buil t 
for use on the 100-m GBT ( | Jewell fc Prestagef [2004] ). 
MUSTANG uses reimaging optics with a pixel spacing 
of 0.63/A, operates in a bandpass of 81-99 GHz, and 
is cooled by a pulse tube cooler and Helium-4/Helium-3 
closed cycle refrigerator. Furt her technical details abou t 
MUSTANG can be found in picker et aD ( [20061 [2008). 
More detailed information is also provided about MUS- 
TANG, the observing strategy, and the data analysis al- 
gorithms in |Cotton et ah] ( |2009l ) and picker et al | j 2009 ) , 
which present the results of other early MUSTANG ob- 
servations. Further information can be obtained at the 
MUSTANG web sitcB 

The observations we present were collected in two runs, 
one on 21 Feb. 2009 and one on 25 Feb. 2009, each ap- 
proximately four hours in duration including time spent 
setting up the receiver and collecting calibration obser- 
vations. For both runs the sky was clear with ~ 6 mm of 
precipitable water vapor, corresponding to ~ 20 K zenith 
atmospheric loading at 3.3 mm. Both were night-time 
sessions, important because during the day variable gra- 
dients in the telescope structure's temperature degrade 
its 90 GHz performance significantly. The telescope fo- 
cus was determined by collecting small maps of a bright 
calibrator source at a range of focus settings; every 30-40 
minutes throughout the session the beam is checked on 
the calibrator. Typically the required focus corrections 
are stable to a few millimeters over several hours once 
residual daytime thermal gradients have decayed. 

Once the focus was established, the in-focus and out-of- 
focus beam maps were used to solve for primary aperture 
wavefront phase errors using the " Out-of-Focus" (OOF ) 
holography technique described by |Nikolic et al.| (|2007). 
The solutions were applied to the active surface. This 
procedure improves the beamshape and increases the 
telescope peak forward gain, typically by ~ 30%. This 
approach is effective at correcting phase errors on scales 
of 20 meters or larger on the dish surface, but is not 
sufficiently sensitive to solve for smaller scales. There- 
fore there are residual uncorrected wavefront errors that 
result in sidelobes out to ~ 40" from the main beam. 
Deep beam maps were collected on the brightest 90 GHz 
sources on several occasions, principally in test runs on 
24/25 March 2009. The repeatability of the GBT 90 
GHz beam after application of the OOF solutions was 
found to be g ood. The analysis of the beam map data is 
discussed in ^ 

(J2000 coordinates 



Maps of 



|3X 

RXJ1347-1145 
13ft47ra30.5s, -11°45'09") were collected with a 
variety of scan patterns designed to simultanously 
maximize on-source time and the speed at which the 
telescope moved when crossing the source. The effects 
of atmospheric and instrumental fluctuations, which 
become larger on longer timescales, are reduced by faster 
scan speeds. The primary mapping strategies were: a) 
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a "daisy" or "spirograph" scan in which the source of 
interest is frequently recrossed; and b) a "billiard ball" 
scan, which moves at an approximately constant speed 
and has more uniform coverage over a square region 
of interest. The nominal region of interest in this case 
was 5' x 5', centered on RXJ1347-1145. The size of the 
maps is sufficiently small that except under the most 
exceptionally stable conditions, instrument and random 
atmosphere drifts dominate the constant atmosphere 
(sec(za)) term and any possible ground pickup. The 
total integration time on source was 3.4 h. 

The asteroid Ceres was observed on both nights and 
used as the primary flux calibrator assuming Tb = 148 K 
(T. Mueller, private comm.). We assign a 15% uncer- 
tainty to this calibration. We checked the Ceres calibra- 
tion on nights when other sources (Saturn, CRL2688) 
were visible and found consistent results to within the 
stated uncertainty. Using these observations and the 
lab-measured receiver optical efficiency of r] op t,rx — 
50 ± 10% we compute an overall aperture efficiency of 
Vaperture,gbt — 20%, corresponding to a Ruze-equivalent 
surface RMS of 315/im. This result is consistent with re- 
cent traditional holographic measurements of the GBT 
surface. Since the observations presented here the sur- 
face has been set based on further holography maps and 
now has a surface RMS, weighted by the MUSTANG 
illumination pattern, of ~ 250/im. 

3. DATA REDUCTION 

3.1. Beam Characterization 

Imaging diffuse, extended structure requires a good 
understanding of the instrument and telescope beam re- 
sponse on the sky. To achieve this we collected numer- 
ous beam maps through our observing runs, including 
several deep beammaps on bright (5 Jy or more) sources. 
After applying the Out-of-Focus holography corrections 
to the aperture the beam results were repeatable; Fig- 
ure [l] shows the radial beam profile from maps of a bright 
source (3C279) collected on two occasions. We find a sig- 
nificant error beam concentrated around the main lobe 
which increases the beam volume from 87 arcsec 2 (for the 
core component only) to 145 arcsec 2 . We attribute this 
beam to residual medium and small scale phase errors 
on the primary aperture. The beam shape and volume 
is taken into account when comparing to model predic- 
tions. By way of comparison, Figure [I] also shows the 
profile of the beam determined from the radio source in 
the center of RXJ1347-1145. Since the SZ map has been 
smoothed, the apparent beam is slightly broader, but al- 
lowing for this, still consistent with the beam determined 
on 3C279. 

3.2. Imaging 

A number of systematic effects must be taken into ac- 
count in the time domain data before forming the image: 

1. The responsivities of individual detectors are mea- 
sured using an internal calibration lamp that is 
pulsed periodically. Optically non-responsive de- 
tectors (10-15 out of 64) are flagged for removal 
from subsequent analysis. Typical detector respon- 
sivities are stable to 2 — 3% over the course of sev- 
eral hours. 
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Fig. 1. MUSTANG+GBT beams determined from observa- 
tions of 3C279. The green dash triple-dot line shows a double- 
Gaussian fit to the observed beam. The purple dashed line is the 
cumulative fractional beam volume. We attribute the excess power 
in the wings of the observed beam to residual medium and small 
scale phase errors on the dish. The data points (diamonds with 
error bars) show a complementary determination of the beam from 
the 5 mJy radio source in the center of RXJ 1347- 1145. The beam 
in this case is slightly wider due to the smoothing (4" FWHM) 
applied to the final map, which is accounted for in the analysis. 

2. Common mode systematic signals are subtracted 
from the data. These are caused by atmospheric 
and instrumental (thermal) fluctuations. The pulse 
tube cooler, which provides the 3K base tempera- 
ture of the receiver, induces a 1.4 Hz signal due 
to small emission fluctuations of the 3K optics. 
The pulse-tube signal is removed by fitting and 
subtracting a 1.41 Hz sine wave. The remaining 
common mode signal is represented by a template 
formed by a weighted average of data from good 
pixels; this template is low-pass filtered and sub- 
tracted from the data, with a fitted amplitude per 
detector. The low-pass filter time constant (typi- 
cally 0.1 Hz) is determined by the stability of the 
data in question. This procedure helps to preserve 
large-scale structure in the maps. 

3. Slow residual per-pixel drifts are removed using 
low-order orthogonal polynomials. 

4. Individual detector weights are computed from the 
residual detector timestreams after the above steps. 
Since the noise level of the detectors varies consid- 
erably this is an important step. Best results are 
obtained by retaining only the top ~ 80% of re- 
sponsive detectors. 

5. The remaining calibrated detector timestreams are 
inspected visually on a per-scan (typically 5 minute 
period) basis. Scans which have timestreams 
with obvious, poorly-removed systematic signals 
remaining are removed. This results in flagging 
28% of scans. The SNR in an individual detector 
timestream is sufficiently low that this does not 
bias our map. 

Following these calibration steps the detector timestream 
data are gridded onto a 2" pixellization in Right Ascen- 
cion and Declination using a cloud-in-cell gridding kernel. 



To check our results we have implemented three, 
mostly independent analysis pipelines. The results in 
this paper are based on a straightforward but flexi- 
ble single-pass pipeline written in IDL, described above. 
There is also an iterative, single-dish CLE AN based ap- 
proach implemented in the OBIT package ( |Cotton[ 2008 ) 
and an optimal SNR method in which the time domain 
data are decomposed into noise (covariance) eigenvec- 
tors; their temporal power spectra computed; and a 
maximum likelihood map constructed from the noise- 
weighted eigenvectors. Results obtained with these al- 
gorithms were consistent. The first two appro aches are 
described in more detail in |Dicker et al.| ( |2009| ). 

Our final map, smoothed by a 4" b WHM Gaussian and 
gridded on 0" .5 pixels, is shown in Figure [2j along with 
the difference between the two individual night maps. It 
shows a strong, clear SZ decrement, well separated from 
the central point source and consistent w ith the level ex- 
pected from the Kitayama et al. ( 2004| hereafter K04) 
150 GHz measurement. The right nand panel shows the 
image formed by differencing the images of the two in- 
dividual nights. By computing the RMS in a fiducial 
region in the center of the difference image (and scaling 
down by a factor of 2 to account for the differencing and 
the shorter integration times) we estimate a map-center 
image noise of ~ 0.3mJy/bm (rms). The noise level in 
regions of the map outside the fiducial region is corrected 
for exposure time variations assuming Gaussian, random 
noise with a white power spectrum. The enhancement 
of the SZE to the south east of the X-ray peak, origi- 
nally detected by Komatsu et al. at 4.2cr significance, 
is confirmed by our measurement at 5.4<j (indicating the 
peak SNR per beam) with a factor of ~ 2 greater angular 
resolution. 

A det ailed assessment of the impact of this is presented 
in § |3.3| Work is underway to develop analysis techniques 
which account for correlated noise in a way that permits 
quantitative model fitting. 

3.3. Simulations 

It is difficult to measure diffuse, extended structure 
such as the SZE, particularly in the presence of po- 
tentially contaminating systematic signals such as time- 
varying atmospheric fluctuations. To assess the impact 
of residual, unmodelled noise fluctuations in the maps we 
have undertaken an extensive suite of simulations which 
replace the raw detector timestream data with simulated 
data. As a source for the simulated data we used real 
detector timestreams collected during observations of a 
blank patch of sky collected for another project. The 
phase of these timestreams with respect to the telescope 
trajectory on RXJ1347-1145 was randomly shifted to cre- 
ate different instances of noisy cluster observations. We 
added simulated astronomical signals as described below 
in order to determine how well the (known) input signals 
are recovered in the maps. 

To assess the spatial fidelity of our reconstructed im- 
ages, random white- noise skies were generated on 2" 
grids, subsequently smoothed by a 9" (FWHM) Gaus- 
sian. These skymaps served as input to generate fake 
timestreams which were then processed by the exact pro- 
cessing scripts used to produce the image in Figure [2] 
The ratio of the absolute magnitude of the Fourier trans- 
form of the reconstructed sky map to the absolute mag- 
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Fig. 2. — MUSTANG image of the SZE in RXJ1347-1145 (left); the same, with the point source subtracted as described in the text 
(center); and the individual nights imaged separately and differenced, on the same color scale (right). The noise in the center of the map 
is is ~ 0.3mJy/bm; contours in the left panel correspond to SNR of 1 to 5 in la increments and account for variations in integration 
time in the map, so are not directly proportional to the image. Color scale units are mJy/bm. The MUSTANG beam (10 /x FWHM after 
smoothing) is shown in the lower right of each panel. 



5 



nitude of the Fourier transform of the input skymap 
measures the fidelity of our image reconstructions as a 
function of angular scale. The results of repeating this 
100 times, with different white-noise skies and noise in- 
stances, are shown in Figure [4j We find that our pipeline 
faithfully recovers structures up to 60", with reasonable 
response but some loss of amplitude on larger scales, up 
to 120". The loss of structure on small angular scales 
is an effect of our relatively coarse pixellization. Simu- 
lations were carried out at similar signal to noise ratios 
as those in our final map, although changes in the signal 
to noise ratio of over a factor of 5 showed no significant 
change to our transfer function 

The common-mode subtraction, essential to remov- 
ing atmospheric and instrumental systematic signals, can 
also introduce negative bowls around bright point sources 
which could mimic the SZE in cases such as RXJ 1347- 
1145. To determine the magnitude of this systematic 
we have followed a similar approach. Instead of white- 
noise skies the input signal consists of a single unresolved 
source with a flux density of 5 mJy at the location of the 
radio source seen in RXJ1347-1145. The resulting nega- 
tive bowl in the reconstructed images has a mean peak 
spurious decrement ~ 2% of the point source peak bright- 
ness, in comparison with ~ 50% for our real data. Ad- 
ditionally the iterative pipeline (OBIT) is much less sus- 
ceptible to such artifacts, and shows consistent results. 
We conclude that this is not a significant contribution to 
our result. 

3.4. Image Domain Noise Estimate 

We divide the data set in half and subtract the individ- 
ual night images to obtain a difference map. The RMS 
of this difference map in the central 85 by 93 arcseconds, 
dividing by two to correct for the differencing and the 
reduced integration time in each individual night image, 
gives an image noise level of 0.3 mJy. A histogram of 
the pixel values in this region of the difference image is 
shown in figure [3j 
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Fig. 3. — Histogram of pixel values in the difference image in Fig- 
ure [2] The histo gram is well described by a gaussian distribution 
of la = 0.3mJy/bm 



We obtain a complementary estimate of the noise from 
a region of the final SZ map well away from the cluster. 



Correcting for the difference in integration times in these 
regions of the map, this result is consistent to within 8%. 

Using this noise figure, the peak SNR per beam in the 
map — on the SZE decrement SE of the cluster core — is 
5.4<j. More aggressive filtering the map results in an even 
higher detection significance for the SE enhancement by 
reducing the low spatial frequency tail of the noise power 
spectrum (Figure |4|. 

3.5. The Effects, and Subtraction, of the Central Radio 

Source 

Our final map has sufficient angular resolution to dis- 
tinguish the central radio source from the structures of 
interest. In particular it is clear that, as seen in ear- 
lier analyse s of the SZE in this cluster ( |Pointecouteau 
et al.|[l999|) there is a strong azimuthal variation in the 
intensity of the SZE at a radius of ~ 20" from the X-ray 
centroid, which also coincides with the radio source. 

To produce a source-subtracted image we fit and sub- 
tract an azimuthally symmetric , do uble- Gaussian beam 
(as determined from 3C279 in § 3.1). The reason for as- 
suming azimuthal symmetry is that the hour angle sam- 
pling of the 3C279 data is considerably more limited than 
that of the RXJ1347-1145 data; therefore the 3C279 data 
will not provide a good measurement of the effective two- 
dimensional beam, only of its average radial profile. Fur- 
thermore the SNR on the point source in RXJ 1347- 1145 
is insufficient to measure significant departures from az- 
imuthal symmetry. 

The average radial profile of the central source in 
RXJ 1347- 1145 is shown in Figure [T]out to r = 15", where 
the signal becomes too weak to measure above thermal 
noise and variations in the SZE. 

3.6. Effect of Background Anisotropics 

The angular scales reconstructed in our map (~ 1'.5 
and smaller) correspond to spherical harmonic multipoles 
of £ = 7200 and higher. On these scales intrinsic CMB 
anisotropics are strongly suppressed by photon diffusive 
damping at the last scattering surface and do not con- 
tribute measurably to our result at the sensitivity level 
we have achieved. 

4. INTERPRETATION & CONCLUSIONS 

4.1. Comparison with Previous SZE Observations 

Figure [5] presents a direct comparison of the MUS- 
TANG and NOBA results in units of main-beam av- 
eraged Compton y parameter. For a more accurate 
comparison, we downgrade the resolution and pixelscale 
of the MUSTANG map to match that of NOBA (13" 
FWHM on a 5" pixel grid). The overall agreement be- 
tween the maps is excellent, in particular as regards the 
amplitude and morphology of the local enhancement of 
the SZE south-east of the cluster core. The largest dis- 
crepancy is south west of the cluster, where NOBA shows 
a 3<j compact decrement which is absent from the MUS- 
TANG data. Considering the low and uniform X-ray 
surface brightness in the vicinity of this discrepancy (see 
Figure [7| and the higher angular resolution and lower 
noise ofthe MUSTANG data, it is likely that this feature 
is a spurious artifact in the NOBA map. Both datasets 
also show a ridge extending north from the shock front 
on the eastern side of the cluster. In the 150 GHz map 
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Fig. 4. — Map noise characteristics. Solid curve: The ratio of the 
absolute magnitude of the Fourier transform of the input map to 
the absolute magnitude of the Fourier transform of the output map 
from simulations. All structure on scales smaller than 1' was recov- 
ered well although there is a fall off towards high spatial frequencies 
due to pixellization effects. Diamonds: The absolute magnitude of 
the Fourier transform of the (reconstructed simulated) cluster map, 
normalized to a peak value of unity. Plus marks: Absolute mag- 
nitude of the fourier transform of signal-free simulated maps, with 
the same normalization as the cluster data. The dashed line shows 
the size of the MUSTANG instantaneous field of view. 
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Fig. 5. — Comparison of NOBA and Mustang maps of RX J1347- 
1145. The color scale shows the MUSTANG data with 5" pixels, 
smoothed to match the published NOBA map resolution. The 
contours show the NOBA map at intervals of 1.16 x 10 -4 in y 
starting from 5 x 10 -5 . Three labeled features are discussed in 
the text: 1, the hot shock south-east of the cluster core; 2, an 
enhancement in integrated pressure to the east; and 3, a compact 
decrement observed at ~ 3cr by NOBA that is absent from the 
MUSTANG image. 

the feature is of marginal significance (1 — 2a); interest- 



Fig. 6. MUSTANG SZE image of RXJ1347-1145 with contours 
(thin green lines) at —1.5, —1.0, and —0.5 mjy/bm. The bold white 
contours, at the same surface brightness levels, show the model SZ 
signal discussed in the text. 



ingly, it is clearly visible in the 350 GHz SZE increment 
map but K04 dismiss it due to the possibility of confusing 
dust emission from the nearby galaxies. 

4.2. Empirical Model of the SZE in RXJ1347-1145 

We construct a simple empirical model for the clus- 
ter SZE assu ming the isothermal /3-model of |Schindler| 
et al.| (|1997[ ) norm alize d by the SZE m e asurem ent of 
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Reese et al.| fl2002| ) and |Kitayama et al | fl2004|) to de- 
scribe the bulk cluster emission. We add a 5 mJy point 
source in the cluster core, coincident with the peak of 
the /3-model, and two Gaussian components in integrated 
pressure, one south-east and one almost directly east of 
the cluster center. In comparing t o our 90 GHz data, 
we use the relatavistic correction of |Sazonov fc Sunyaev 
( fl998| ), assuming kT = 25keV (which reduces the am- 
plitude of the decrement by 15%) for the Gaussian com- 
ponents and kT = lOkeV for the bulk component. The 
parameters chosen (two Gaussian widths for each com- 
ponent, a position, a peak surface brightness, and a posi- 
tion angle) are shown in Table IT] T he resulting sky image 
is conv olved with our PSF (§ K3Jl) and transfer function 
(§ 3.3). We find that this provides a good match to the 
data (Figure [6]). The peak comptonization at 10" Gaus- 
sian resolution is 3.9 x 10 -4 on the eastern ridge and 
6.0 x 10 -4 on the region identified as a shock by Ko- 
matsu et al. When convolved to 19" FWHM (NOBA) 
resolution, we find Ay — 3.9 x 10 -4 , close to their ob- 
served value Ay — 4.1 x 10 -4 . The intent of this static, 
phenomenological model is simply to provide a descrip- 
tion of the observed high angular-resolution SZE and 
a direct comparison of NOBA and MUSTANG results. 
Work is underway which will allow quantitatively deter- 
mining the best fit physical model by simulataneously 
fitting datasets at multiple wavelengths using a Monte- 
Carlo Markov Chain. This work is beyond the scope of 
this paper and will be presented in a follow up publica- 
tion. 
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Fig. 7. — Left: False color composite Image of RXJ1347-1145. Red/blue: Mustang SZ. Green: Archival HST/ACS image tak en through 
the F814W filter; and white contours: Surface mass density k from the weak + strong lensing analysis of |Bradac et al.| ( |2008| >. Contours 
are linearly spaced in intervals of Ak = 0.1 beginning at k = 1.0. Several features are labelled: A indicates the central BUG, which is 
a radio source; B indicates the BCG of the secondary cluster; and 3 (also labeled in the right-hand panel) indicates the location of the 
discrepancy between NOBA and MUSTANG, discussed in the text and Figure [5] The diamond, cross and box mark the locations of the 
peaks in X-ray surface brightness, surface mass density and SZE decrement respectively. Right: Contours of the MUSTANG decrement 
SNR (la to 5cr in la increments) superposed on the Chandra count-rate image smoothed to 10 /x resolution. 

Component Amplitude Offset Notes 

[y/io- 3 ] n 

/3-model 1.0 0,0 C = 10", (3 = 0.60 

Shock 1.6 -14, 14 ai = 8", a 2 = 2", P.A. = 45° 

Ridge L0 10, 14 ai = 8", a 2 = 2", P.A. = -15° 



TABLE 1 

Note: Offset is (north, east) of peak X-ray position 



4.3. Multi-wavelength Phenomenology 

Our data show an SZE decrement with an overall sig- 
nificance of 5.4cr. At the center of the cluster, coincident 
with the peak of X-ray emission and the brightest cluster 
galaxy (BCG), there is an unresolved 5 mJy radio source. 
This flux density is consistent with the 90 GH z flux den- 



sity presented in |Pointecouteau et al. (2001), as well as 
what is expected from a power law extrapolation of 1.4 



GHz and 30 GHz measurements dCondon et al. 1998 



Reese et al. 2002). A strong, localized SZE decrement 



can be seen 20" to the south-east of the center of X-ray 
emission and clearly separated from the cluster center. 
Our data also indicate a high-pressure ridge immediately 
to the east of the cluster center. 

K04 tentatively attribute the south-east enhancement 
to a substructure of gas 240 ± 183 kpc in length along 
the line of sight, at a density (assumed uniform) of 
(1.4 ± 0.59) x 10 -2 cm -3 and with a temperature — 



28.5=b7.3ke V. Recent X-ray spectral measurements ( |Ota 
|et al. ||2008l ) with SUZAKU also indicate the presenceof 
hot gas in the south-east region (T e = 25.I+45 +g'J keV 
with statistical and s ystematic er r ors, r espectively, at 
90% confidence level). | Allen et al. (2002) have reported 
that the slight enhancement ol sorter X-ray emission in 
this region seen by Chandra is consistent with the pres- 
ence of a small substructure of hot, shocked gas. |Ki- 



|tayama et al.| ( |2004[ ) attribute the hot gas to an ongoing 
merger in the plane of the sky. The merger hypothe- 
sis is supported by optical data, in particular, the pres- 
ence of a second massive elliptical ~ 20" directly to the 
east of the BCG that coincides with the center of X- 
ray emission (and with the radio point source). Further- 
more the density and temperature of the hot substruc- 
ture indicate that it is substantially overpressured com- 
pared to the surrounding ICM. Assuming a sound speed 
of 1600 km/sec this overpressure region should relax into 
the surrounding ICM on a timescale ~ 0.1 Gyr, again 
arguing for an ongoing merger. 

Our data support this merger scenario. To put them in 
context, Figure [7] shows a composite image with archival 
Chandra and HST data, and the weak + strong lensing 
mass map of Bradac et al. (2008). We propose that the 
data are best explained by a merger occuring in or near 
the plane of the plane of the sky. The left-hand ("B") 
cluster, having fallen in from the south- west, has just 
passed closest approach and is hooking around to the 
north-west. As the clusters merge shock forms, heating 
the gas in the wake of its passage. As argue d by |Ki- 
taya ma et al. (2004), and seen in simulations (Takizawa 
1999), the clusters must have masses within a factor ol 
2 or 3 of equality and a substantial (~ 4000 km/sec) 
relative velocity in order to produce the high observed 
plasma temperatures, T e > 20 keV. This merger geome- 
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try is consistent with the lack of structu re in the line-of- 
sight cluster member galaxies' velocities (Cohen & Kneib 
[2002]) . The primary (right-hand, "A") cluster contains 
significant cold and cooling gas in its core (a "cooling 
flow" ) . Such gas is see n to be quite robust in simulated 



majo r cluster mergers ( Gomez et al 
Even in cases where the coo' 



2002 



mg 



Poole et al 
n ow is finally 



2008}. 

disrupted by the encounter, | Gomez et al.| ( |2002[ ) find a 
delay of 1 — 2 GYr between the initial core encounter 
and the collapse of the cooling flow. The existence of 
a strong cooling flow, therefore, does not argue against 
a major merger in this case. More detailed simulations 
could shed further light on this interesting system. 

4.4. Broad Implications 

Since calibrating SZ observable - mass relationships is 
vital to understand the implications of ongoing SZE sur- 
veys, it is important to understand the mechanism by 
which a substantial portion of the ICM can be heated so 
dramatically, and how this energy is distributed through 
the ICM ove r time. Observations o f cold fronts in other 
clusters (e.g. Vikhlinin et al.p OOl) have shown that en- 
ergy transport processes in the ICM are substantially in- 
hibited, perhaps by magnetic fields. It is distinctly pos- 
sible, then, that once heated by shocks, very hot phases 
would persist. 

We have estimated the magnitude of the bias in an 
arcminute-resolution Compton y measurement that is in- 
troduced by a hot gas phase by convolving the two gaus- 
sian components of the SZE model in Table [T] with a V 
FWHM Gauss ian beam, typical of SZ survey telesco pes 
such as ACT ([Swetz et al.|[2008l |Hincks et al. 200 9]) o r 
SPT jstaniszewski et al ||20lJgf[Vanderlmde et al ||2U10[ ). 
Compared with the bulk emission component, also con- 
volved with a V beam, the small-scale features are a 10% 
effect. While relatively modest this is a systematic bias 
in the Compton y parameter which, if not properly ac- 
counted for, would result in a 20% overestimate in dis- 
tances (underestimate in Ho) derived from a comparison 
of the SZE and X-ray data which did not allow for the 
presence of the hot gas component. To assess the impact 



on the scatter in M — y, a larger sample of high-resolution 
SZE measurements is needed. A full calculation would 
also need to take into consideration effects such as de- 
tection apertures and the spatial filtering due to imaging 
algorithms, some of which would increase the importance 
of the effect and some of which would decrease its impor- 
tance. 

This is the one of a very few clusters that has been ob- 
served at sub-arcm inute resolution in the SZE (see also 



|Nord et al. 
exhibit simi 



2009), so it is possible that many clusters 
ar behavior. 



Such events, if their enhance- 
ment of the SZE brightness is transient, could also bias 
surveys towards detecting kinematically disturbed sys- 
tems near the survey detection limits. The astrophysics 
that has been revealed by high resolution X-ray obser- 
vations, and is beginning to be revealed by high resolu- 
tion SZE data, is interesting in its own right. The SZE 
observations require large-aperture millimeter telescopes 
which have henceforth been lacking, but with both large 
single dishes and ALMA coming online, exciting observa- 
tions will be forthcoming. There is substantial room for 
improvement: since the observations we report here the 
GBT surface has improved from 320/im RMS to 250/im 
RMS, which will yield more than a factor of 1.5 improve- 
ment in sensitivity. The array used in these observations, 
while state of the art, has not yet achieved sky photon 
noise limited performance; further progress is being made 
in this direction. Considering these facts, and that the 
results presented here were acquired in a short period of 
allocated telescope time (8h), this new high-resolution 
probe of the ICM has a bright future. 

The National Radio Astronomy Observatory is a facil- 
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map; Masao Sako, Ming Sun, Maxim Markevitch, Tony 
Mroczkowski and Erik Reese for helpful discussions; and 
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on the manuscript. 
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